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Abstract. The dissociative recombination of the lowest rotational states of has been 
investigated at the storage ring TSR using a cryogenic 22-pole radiofrequency ion trap as 
injector. The was cooled with buffer gas at ~15K to the lowest rotational levels, 
(J, G) = (1,0) and (1,1), which belong to the ortho and para proton-spin symmetry, respectively. 
The rate coefficients and dissociation dynamics of H^{J, G) populations produced with normal- 
and para-iif2 were measured and compared to the rate and dynamics of a hot beam from a 
Penning source. The production of cold rotational populations was separately studied by 
rovibrational laser spectroscopy using chemical probing with argon around 55 K. First results 
indicate a ~20% relative increase of the para contribution when using para-iif2 as parent gas. 
The rate coefficient observed for the para-i/2 source gas, however, is quite similar to the 
rate for the normal-//2 source gas. The recombination dynamics confirm that for both 
source gases, only small populations of rotationally excited levels are present. The distribution 
of 3-body fragmentation geometries displays a broad part of various triangular shapes with 
an enhancement of ~12% for events with symmetric near- linear configurations. No large 
dependences on internal state or collision energy are found. 



1. Introduction 

The ion is of utmost theoretical and astrophysical importance. It is the simplest polyatomic 
ion and thereby constitutes an important benchmark ion for the theory of both molecular 
structure and molecular collision dynamics. It is also a key ion in the universe as it initiates 
a network of reactions leading to the formation of complex molecules. The dissociative 
recombination (DR) of plays a crucial role as destruction mechanism for the ions 
[U [2] and as production mechanism for the neutral hydrogen fragments. The DR reaction 
has been extensively researched over the past decades [3l S] and with the knowledge on the 
reaction increasing, the focus is moving to more and more detailed investigations. Internal-state 
dependences of the DR reaction, at zero as well as at elevated electron energies, are now being 
studied O E]. Disagreements on the magnitude of the DR rate coefficient of have been 
dominating the scene for quite some time 0. Presently, the differences that have been observed 
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are attributed to different internal-state distributions of the ions and consent between storage- 
ring experiments has been reached [U [9]. Furthermore, previous discrepencies with afterglow 
studies seem to have been resolved recently and the latest theoretical calculations show good 
overall agreement with experiment |10 t lll t [T2]. However, the predicted rate coefficient still shows 
local discrepancies up to an order of magnitude around several electron energy regions. Again, 
differing internal-state distributions are given as possible explanation. In the cold interstellar 
space, the difference of the DR rate between the two lowest rotational levels, (J, G)=(1,0) and 
(1,1) at an energy difference corresponding to kT=33K only, may be of particular importance. 
Since for these two levels are equivalent with the ortho and the para symmetry of the total 
nuclear spin, respectively, we shortly denote the two levels as ortho- and para-if^. Previous 
experiments indicate that the DR rate coefficient of para-H^ could be larger than that of ortho- 
for specifically the low electron energies prevailing in the cold interstellar medium [9j. In 
these experiments, para--ff2 was used as precursor gas and although the rotational populations 
of the ion beam could not be measured, an enhancement of the para--?/^ population is 
expected in this case. The higher para-i/^ rate was until recently in contradiction to theory, 
which predicted the opposite [HI [13]. However, new theoretical results now state a roughly ten 
times higher rate for para-H^ at low temperature |12j . 

The DR of cold has two fragmentation channels for low-energy electrons, a 3-body 
breakup containing only ground-state atoms and a 2-body breakup that includes the H2 product 
molecule, 

H+{J, G) + e-{Ed) ^ H{ls) + H{ls) + H{ls) , 

^ H2{v,J)+H{ls) 

Here, the cold initial ion is given to be in a low rotational state of its ground vibrational state 
with J the total angular momentum associated with the motion of the nuclei and G = \K — l\, 
where K is the projection of J onto the molecular axis and I is the vibrational angular 
momentum. The product atoms are in their ground state, whereas the H2 molecule may be 
rovibrationally, (f, J), excited. Further dissociation channels become possible when additional 
energy, Ed, is introduced via the electrons. At typical molecular-cloud temperatures of 10-60 K, 
populates the two to three lowest rotational levels, (J, G)=(1,0), (1,1), and (2,2), of the 
vibrational ground state, where the (2,2) level belongs to the para spin symmetry. 

We present a two- fold investigation of cold . The first focus is on the production and 
characterisation of rotationally cold populations using normal--fr2 and para-if2 as parent 
gas. In order to probe the different rotational distributions, rovibrational spectroscopy on the 
ions in a 22-pole radiofrequency ion trap is implemented. The second focus is on the DR 
of H^{J,G) populations with altered ortho-to-para ratio. High-resolution rate coefficients up 
to 0.4 eV are given and additionally compared to a hot beam to investigate rotational 
heating. The issue of internal excitation is addressed by fragment imaging. For the first time, 
the dissociation of were studied using the 3D imaging technique [13], and 3-body breakup 
from both 2D and 3D imaging are presented. 

2. Experimental Setup 

The cold samples are prepared in a 22-pole radiofrequency (rf) ion trap setup, as shown in 
the inset of Figure [H The ions are produced in a rf storage ion source [15] through electron- 
impact ionisation of H2 and subsequent collisions of the H2 with the hence produced H2 ■ The 
ions are then guided through a rf quadrupole to the 22-pole rf ion trap [16] . Here, the ions 
are trapped radially by the 22 poles and axially by an entrance and exit electrode. The number 
of stored ions may be varied from a few to several millions of ions, whereas the storage time may 
be varied from ms to several seconds. During storage, the ions are buffer-gas cooled with helium 



at a chosen ambient temperature of 10-60 K. For more details on the 22-pole trap setup, the 
reader is referred to Refs. [l7l[l8]. In the present work, the setup is used as a stand-alone device 
for probing the H^{J,G) populations implementing rovibrational spectroscopy and it is used 
as an ion injector for the storage ring TSR to the study of the DR of the cold (see Figure 
[l]). Unfortunately, the conditions required for the spectroscopy measurements are different from 
those of the DR measurements. As such, care has to be taken when deploying information from 
the laser-probed populations for the TSR measurements. 




Figure 1. The 22-pole trap setup used as injector for the ion storage ring TSR. The ions 
are prepared in the trap, accelerated to 4.05 MeV, and injected into the ring. A cold electron 
beam, produced with a liquid-nitrogen cooled photocathode [TH], is merged with the ion beam 
in the electron target. The neutral fragments are detected either with a surface barrier counting 
detector or a position-sensitive imaging system [2^ I21j. 



2.1. Spectroscopy Experiments 

To verify the buffer-gas cooling inside the 22-pole trap and to probe the rotational-state 
populations, an action spectroscopy scheme is used. To this purpose, rovibrational spectroscopy 
in the infrared is implemented as lacks any stable electronically excited states. As a sensitive 
probe of the vibrational transitions, a laser-induced chemical reaction with argon as reactant is 
used. The charge exchange reaction, 

H^{v = 0) + Ar ^ ArH+ + H2- 0.57 eV (1) 

is endothermic by 0.57 eV for ground-state ions. However, if at least two vibrational quanta 
are excited in , the reaction becomes energetically possible. Once the ions in the trap are 
cooled down to the chosen temperature, a tunable infrared diode laser of ~1400 nm is used to 
trigger transitions from the lowest rotational states of to the third harmonic of the bending- 
mode vibration. The laser-excited ions will react with the abundant argon in the trap to form 
ArH^ , which is trapped as well and used as spectroscopic signal. The ions are extracted from 
the trap, analysed by a mass spectrometer, and directed to a scintillation detection system 
capable of counting single ions with near-unity efficiency. These spectroscopy measurements are 
performed at a minimum temperature of 55 K to keep the argon in the gas phase. For more 
details, see Ref. [T7j . 

The action spectroscopy scheme has been used before to probe cold produced with 
normal- i?2 [iZ]- Here, also the ^^3" rotational distribution produced with para-i?2 is presented. 



Moreover, important improvements have currently been implemented to enhance the cooling of 
the in the trap as well as to enhance the ArH~^ spectroscopy signal. Gas purities, laser 
power (increased by a factor of 10), and frequency resolution have now been refined. As a 
result, the typical ArH~^ count at resonance is 1.5-2 per trap filling of 300 ions, an order 
of magnitude higher than before. The number densities of the helium and argon gas inside the 
trap have been optimised to ~10^^ cm~^ and ~10^^ cm~^, respectively. The hydrogen number 
density, ~10^'' cm~^, and the trap storage time, 250 ms, were unchanged. 

2.2. Dissociative Recombination Experiments 

The DR measurements were performed with the 22-pole trap as injector of the cold ion beam 
for the storage ring TSR. The trap temperature was set to 15 K in order to simplify the DR 
measurements, as only the lowest para and ortho rotational states, (1,0) and (1,1), are populated 
at that temperature. The number of ^^3" ions had to be hugely increased to the maximum trap 
filling of around 2 • 10® in order to supply a sufficient DR count rate [22]. The ortho-to- 
para ratio is dominantly determined by collisions with H2, therefore the H2 number density 
(10^*^ cm~^) was kept the same as in the spectroscopy experiments. The helium number density 
was increased to ~10^^ cm~^ to compensate for the slightly reduced trap storage time of 100 ms. 
After this cooling time, the ions were released to be accelerated and stored in the TSR at an ion- 
beam energy of 4.05 MeV. The transmission of the ion beam was about 10%, resulting in very low 
counting rates (in a later storage-ring experiment the transmission could be improved to 50%). 
The ring storage time was set to 10 s during which the ion beam was merged in the electron target 
with an ultracold electron beam produced by a liquid-nitrogen cooled photocathode p9]. The 
transversal and longitudinal temperatures of this electron beam were T_l = 1 meV and Ty =30 
^eV, respectively. The electron density was roughly 9 • 10^ cm~'^. Rate-coefficient curves were 
measured with a new large surface barrier detector of 10x10 cm^ at collision energies between 
and 0.4 eV by varying the electron detuning energy, E^. The dynamics of the DR reactions 
were investigated at E^ = and 6 meV through the 2D [20] and, as a first for , the 3D 
|14i I21| imaging technique. The electron cooler was set to eV collision energy to continuously 
cool the ion beam and, for the first 0.5 s of storage, the electron target was additionally set 
to eV in order to speed up the initial translational phase-space cooling of the beam. This 
initial cooling time was deliberately kept low in order to minimise possible internal heating- 
effects from the toroidal regions, where the electron beam is no longer parallel to the ion beam 
and elevated detuning energies may cause rotational excitation of the ions. Only for the 0-eV 
imaging measurements, the electron cooler was switched off after this time, whilst the target 
remained at eV. 

3. Results 

3.1. Rovibrational Spectroscopy 

Figure [2] shows the observed ArH'^ signal for the transitions from the three lowest rotational 
states present at 55 K to the third vibrational level. The upper figure results from ions 
produced with normal- 5 the lower from ions produced with para-ff2- We are currently 
still studying possible effects of the change of the precursor gas on the absolute normalisation of 
the signal. The relative normalised intensities for a given normal- -ff2 or paTa-H2 measurement, 
however, are reproducible and well under control. The intensities are normalised to the laser 
power of ~12 mW, to a trap filling of 300 ions, and to an ArH~^ lifetime of 25 ms, such that 
they may be coarsely compared to the previously published spectroscopy measurement [T7] . 

The experimental changes mentioned before have now lead to improved signal strengths 
giving rise to nicely observable Gaussian profiles. The cooling of the ions also proved to be 
more efficient at the optimised helium pressure and much lower temperatures are observed than 
in 2004 [T7]. The profile widths of the respective transitions denoted by R{J,G) - where R 



Figure 2. The measured absorp- 
tion profiles for the three observed 
transitions in the (0, 3^) <— (0, 0'') 
vibrational overtone band of ^^3" 
produced with normal- -fr2 (upper 
panel) and para-i72 (lower panel) . 
The signal is derived by averaging 
over many trap fillings. 

stands for the (J+l) branch and J,G stand for the rotational quantum numbers of the lower 
state - are determined by Doppler broadening, which conforms a temperature of ~67K for all 
transitions. The line- intensity ratios are determined by the given laser intensity, the Einstein 
coefficients, and the initial-level populations. Assuming the theoretical Einstein coefficients of 
8.9, 4.9, and 8.2 s~^ for R{1, 0), i?(l, 1), and R{2, 2), respectively [23], the populations and hence 
ortho-to-para and para-to-para ratios can be extracted. The "para rotational temperature" is 
derived from the ratio between the two para transitions, R{1, 1) and i?(2,2), and matches the 
translational temperature of ~67K for both ion-source gases. Thermal equilibrium seems to 
have been reached in both cases. The "ortho-to-para temperature" is derived from the ratio 
between the ortho transition, i?(l,0), and the para transitions, i?(l, l)+i?(2, 2). A significant 
decrease of the ortho- -ff^ line relative to the para lines is visible when using para--ff2 as parent 
gas. The ortho-to-para temperature correspondingly decreases, namely from ~140K to ~36K. 
This clearly shows that under certain circumstances the ortho-to-para ratio can become non- 
thermal and can deviate significantly from the translational and rotational temperature of an 
ion sample. The population for the normal-//2 parent gas conforms an ortho-to-para ratio of 
60:40 using the above mentioned Einstein coefficients, whereas the sample produced with 
para-i/2 shows the opposite, an ortho-to-para ratio of 40:60. It seems that the use of para-i?2 
as parent gas leads to enrichment of the rotational states belonging to the para-iif^ symmetry. 

3.2. Dissociative Recombination Rate Coefficients 

As mentioned before, the 22-pole trap was used as an injector for the storage ring TSR, where 
the DR rate of the cold populations was investigated between and 0.4 eV collision energy. 
Also, to investigate rotational heating, the DR rate of a hot beam produced with a Penning 
source was measured. Due to low count rates for the cold ion beams (only 10% transmission), 
insufficient statistics could be gathered around the 10-eV resonance. Therefore, the rate curves 
have been normalised to the 12-meV resonance. This resonance is at an energy low enough for 
sufficient statistics, however, presumably at an energy high enough for the normal- and para-H^ 
rates to be the same, as observed in Ref. [9]. The hot rate curve, which does contain enough 
statistics, is normalised to the 10-eV resonance. Both normalisations are to the photocathode 
measurement of Ref. [9j as it bears a similar energy resolution. Figures [3] and S] show the reduced 
rate coefficients, a ■ with a the energy-dependent rate coefficient that is proportional to 
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^fEd- Any structures observable on the reduced rate are related to resonant behaviour. Figure 
[3] (upper frame) displays the measured rates of the hot and cold H.'^ produced with normal- 
Overall, the hot H'^ rate is higher than the cold K'^ rate, except for the resonance around 7 
and 270 meV. Although the temperature of the H'^ beam from the Penning source is on the 
order of 10^ K (see Section FS-Sp . resonant behaviour is still clearly visible. Figure [3] (lower frame) 
shows the theoretical predictions for 100 K and 1000 K as calculated by Ref. [12]. Similar to 
the experiment, the hot K'^ rate is overall higher, however, at very low energies, the opposite 
behaviour is predicted. Several structures appear both in theory and in experiment. The 
intensities of the theoretical rates seem to be systematically lower than those measured, leading 
to a better agreement between our cold beam and the 1000 K theoretical rate. This raises 
the question whether heating effects in the storage ring might have considerably warmed up the 
initially cold ion beam. The fragmentation dynamics presented in the next section, however, 
demonstrate that this is not the case. 




Figure 3. The upper frame shows the 
measured reduced rate coefficients of the 
cold (black) and hot (grey) produced 
with normal-ff2 in the 22-pole trap and the 
Penning source, respectively. The lower frame 
shows the reduced theoretical rates for 100 K 
(black) and 1000 K (grey), folded with the 
experimental electron temperature from Ref. 



Figure 4. The upper frame shows the high- 
resolution reduced rates of H'^ produced with 
para-ff2) present data (grey), and normal- 
i?2! present data (black) and previous data 
(light grey) 0. The lower frame shows the 
previously measured lower-resolution reduced 
rates of produced with para-i/2 (grey) and 
normal- i72 (black) [S]. 



Figure H] displays the reduced rates of the cold ion beams produced with normal- and 
para-i?2 (upper frame). The H'^ rate using para--ff2 exhibits slightly more structure below 5 
meV collision energies. Otherwise, the rate coefficients are quite similar and show no sign of 
the predicted order of magnitude difference in DR rate between para--ff^ and ortho--ff^ [12]. 
The two populations in the storage ring produced with normal- and para--ff2 are unknown 
and may very well differ less than was the case in the spectroscopy measurements, however, a 
DR rate that is an order of magnitude faster for para-iJ^ than for orXhio-K'^ should manifest 
itself even at small increases in the ^w^-H^ contribution. No large dissimilarities between the 
two source gases are observed that would point towards the predicted ten times more rapid 
dissociation of para-i?^. Still, nearly identical H^{J,G) populations in the storage ring for 
both source gases cannot be excluded. 

Also shown in figure H] is the previous photocathode measurement for normal- ^^3" (upper 




Figure 5. The distributions of the beams for (a) the cold produced with normal- 
and para-if2 in the 22-pole trap and (b) the hot produced with normal- in the Penning 
source. All black curves are Monte-Carlo simulations of isotropic dissociations at the given 
temperatures, considering the allowed rotational levels and their probability in a Boltzmann 
distribution. All distributions are normalised to unit area. The deviation of the measured 
distributions of produced with normal- from the simulations are plotted in the respective 
lower frames. 



frame) and the previously measured rates for the normal- and para-i72 source gases 
as measured with the thermal cathode (lower frame) [H]. The photocathode data clearly 
exhibit more pronounced resonances as the resolution of the liquid-nitrogen cooled photocathode 
(Tg _|_ = 1 meV) |19] is much better than that of the thermal cathode (Tg^^ = 4 meV). A small 
resonance is revealed around 17 meV and below 5 meV structure is unfolded. The increased 
structure visible for the para-i72 source gas raises the rate slightly, in accordance with the 
previous thermal-cathode experiment, i.e., assuming the normalisation on the 12-meV resonance 
is valid. Further analysis is still required. It is noted that the internal distributions in the ring 
may also differ between the present and previous DR experiments, both in ortho-to-para ratio 
for the produced with para-i^2 as in the achieved internal temperature that may give rise 
to a differing contribution from excited rotational levels irrespective of source gas. 

3.3. Dissociative Recombination Dynamics 

The breakup dynamics of the DR of and their possible (J, G) state-dependences have been 
studied as well. Imaging data have been taken for both the 2-body and 3-body fragmentation 
channel. Here, only dissociation into three H atoms is treated. Figure [5^) shows the projected 
distance distributions of our cold beams from the 22-pole trap as acquired by the 2D 
imaging [20] , normalised to unit area. The distance plotted is R'^ = R^ + R2 + -R3 with i?j the 
distance from particle i to the centre-of-mass in the 2D detection plane. Forward Monte-Carlo 
simulations for isotropic dissociations at three different rotational temperatures, normalised 
to unit area, are plotted to reveal a rough estimate of the temperature of the ion beam and 
anisotropy effects. To clarify differences, the deviation of the distribution for the normal-i^2 
source gas from the simulations is displayed in the lower frame. The discrepancy at the peak 
of the distribution could indicate anisotropy and further simulations are needed to investigate 
the observed feature. Figure Eb) shows the R^ distribution of the hot from the Penning 
source together with forward simulations at three high temperatures. Below, the difference of 
the measured data with the simulations is shown. Again, no large anisotropics can be observed. 

The best temperature fit is determined from the least-squares minimum between the falling 
slope of the data and the simulation, with the simulation normalised to the maximum of 




Figure 6. The 3-D imaging data of hot from the Penning source for a) aU dissociations and 
b) the near-transversal dissociations, i.e., dissociations with cos{9)>0.8, where 9 is the angle 
between the normal of the molecular plane and the beam axis. 



the measured data (not displayed) in order to disregard the anisotropy effects. The hot 
distribution appears to have a temperature of around 3000 K. For the cold H^, the least-squares 
deviation is smallest at the lowest temperature and increases with rising temperature. Although 
a significant contribution of the J=2 level cannot be excluded, the temperature is most certainly 
below room temperature, which demonstrates that the ion beams are indeed cold. The apparent 
difference between experiment and theory at low temperatures can therefore not be explained 
by considerable heating of the ion beam. 

Figure [6^) gives the Dalitz plot of the hot beam from the Penning source for 0-eV 
collision energy, for the first time acquired by the 3-D imaging technique [Ti] . The geometries of 
the corresponding dissociations are displayed in the lower sextant. As proof-of-principle for the 
3D imaging, the events dissociating near-perpendicular to the detector are displayed in figure 
Eb)- Near-perpendicular events with cos{9)>0.8 are selected, where 9 is the angle between the 
normal of the molecular plane and the beam axis. The tranversal 3-D Dalitz plot is indeed 
similar to the total 3-D Dalitz plot. The 3-D imaging data are also consistent with the Dalitz 
plots obtained previously by Monte-Carlo reconstructions of 2D projected distance distributions 

m- 

To investigate the distributions for the different populations in more detail, the Dalitz 
coordinates have been transformed to polar coordinates (as illustrated in figure [6^) to be 
able to integrate over the angle and the radius. The integrated radius describes the amount 
of (non-)linearity and the integrated angle the deformation from the isosceles/equilateral 
triangle. Figures [7^) and b) display the integrated-coordinate distributions at eV. A uniform 
background, i.e., isotropic distributed geometries, can be observed superimposed with ~12% 
anisotropic dissociations, which exhibit a preference towards linear breakup, independent of the 
population. As discernable in the Dalitz plot, the near-linear dissociations peak around the 
symmetric linear breakup and drop to (even slightly below) the level of the isotropic background 
for the asymmetric linear geometries. The integrated-coordinate distributions show no significant 
internal-state dependence. Figures [7b) and d) display the distributions for the hot and the cold 
-^3' at 6 meV together with the 0-eV distributions of the hot for comparison. The breakup of 



the cold produced with para--?/2 at 6 meV was not measured. The fragmentation dynamics 
of the hot at 6 meV resemble those at eV. The cold breakup exhibits a small difference, 
peaking slightly more around the symmetric linear dissociations. No large energy dependence 
is observed. 
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Figure 7. The projected Dalitz polar coordinates at a,b) eV and c,d) 6 meV for the given 
populations. For comparison, the hot data at eV are displayed again with the 6 
meV distributions. As illustrated in figure [H ranges from the equilateral triangle to linear 
dissociation geometries and (po ranges from asymmetric to symmetric geometries. To improve 
clarity, only the errorbars of the distribution with the lowest statistics are shown. 



4. Conclusions 

The rovibrational spectroscopy shows that the use of para-/f2 as parent gas can indeed enrich 
the para contribution of the population. In the storage-ring experiment, where the para 
enrichment cannot be confirmed, no increase in the DR rate has been found when using para- 
H2 as source gas that may support the prediction of a 10 times higher paica-H^ rate than 
ortho-ff^ rate. The trend of increasing overall rate with increasing temperature as predicted 
by theory is, however, confirmed by the present data. Furthermore, the ions prepared 
in the 22-pole trap and transferred to the TSR are really rovibrationally cold, featuring a 
temperature that is distinctly lower than 300 K. As such, heating effects cannot account for the 
discrepancy in intensity between the theoretical and experimental rate-coefficient curves. The 
3-body fragmentation dynamics of the DR of show no large dependences on the internal- 
state distribution of the initial ^^3^ nor on electron energy. A fraction of 12% of the ions are 



found to dissociate in symmetric near-linear geometries, while 88% have a uniform distribution 
of geometries. 
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